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Pion and kaon physics with improved staggered quarks* 
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We compute pseudoscalar meson masses and decay constants using staggered quarks on lattices with three 
flavors of sea quarks and lattice spacings «0.12 fm and ~0.09 fm. We fit partially quenched results to "staggered 
chiral perturbation theory" formulae, thereby taking into account the effects of taste-symmetry violations. Chiral 
logarithms are observed. From the fits we calculate and /x, extract Gasser-Leutwyler parameters of the chiral 
Lagrangian, and (modulo rather large perturbative errors) find the light and strange quark masses. 



The masses and decay constants of light pseu- 
doscalars can be computed with high precision at 
fixed lattice spacing and quark mass. Assuming 
that the chiral and continuum extrapolations can 
be performed with controlled systematic errors, 
such a calculation therefore can provide an accu- 
rate determination of the physical properties of 
the tt-K system, as well as a sensitive check on 
our algorithms and methods. 

Our simulations have three light dynamical fla- 
vors with the Asqtad actionpp, with 0(a s a 2 ) 
errors in general, but 0(a 2 a 2 ) taste- violations, 
which require 2-gluon exchange. One group of 
lattices ("coarse") has a ~ 0.12 fm, the dynamical 
strange quark mass held fixed at a nominal value 
am s = 0.05, and five different values of dynami- 
cal light quark mass m u = rrid = m q . The second 
group ("fine") has a » 0.09 fm, a nominal value 
am s = 0.031, and two values of m q . The nom- 
inal strange quark mass turned out to be about 
18% (10%) larger than the physical mass, ™P hys , 
on the coarse (fine) lattices. Table ^ shows the 
lattice parameters used. 

The relative scale among coarse (or, separately, 
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6.81 


20 3 x 64 


262 
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6.76 


20 3 x 64 


447 
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7.09 
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Table 1 

Lattice parameters. The lattice sets above the 
double line are "coarse;" those below are "fine." 

fine) lattices is kept fixed using the length r% |2I3| 
from the static quark potential. We reduce statis- 
tical fluctuations in r\ja by using a "smoothed" 
r\ja coming from a 4-parameter fit to a smooth 



function of 10/g and mtot = 2m q 
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absolute scale is taken from T 2S'-1S' or IP- IS* 
splittings, determined by the HPQCD group [H] 
on the lattices with m q = 0.2m s . We extrapo- 
late to the continuum linearly in a s a 2 and get 
r*i = 0.317(7) fm. Since the physical n has some 
(small) dependence on sea quark masses, fixing 
a using r\/a can introduce some spurious depen- 
dence on mass 0]. This has a negligible effect 
on quantities like /„., but could be a significant 
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systematic for the Gasser-Leutwyler parameters 
(Li). Estimates of this effect are in progress. 

On the coarse (fine) lattices, quark propagators 
are computed every 6 time units for 9 (8) values 
of valence mass between 0.1m s (0.14m s ) and m s . 
We determine meson masses and decay constants 
from simultaneous hts to wall-point and point- 
point correlators. The complete covariance ma- 
trix of decay constants and masses on each lattice 
set is computed using single elimination jackknife; 
we then adjust for autocorrelations by scaling the 
matrix by a factor estimated by blocking lattices 
in groups of four. 

We hrst fit squared masses of various-taste pi- 
ons to the tree-level form [H] : linear in quark mass 
plus a constant splitting (at fixed a) for each taste 
multiplet. Although the fits are poor, they do 
give the masses within ~ 5%. The slopes and 
splittings determined may then be consistently 
used as inputs to the one-loop terms in the chiral 
log fits below. The ratio of splittings on the fine 
lattices to those on the coarse lattices is » 0.35. 
This is consistent with expectation that taste vi- 
olations go like 0(a 2 a 2 ): using a s = ay(q*) at 
one-loop [7] and q* = 3.33/a, this ratio is 0.375. 

We then compare the (partially quenched) de- 
cay constants and squared meson masses to one- 
loop chiral log forms (including finite volume cor- 
rections) computed in staggered chiral perturba- 
tion theory (SXPT) 0. Up through NLO, we 
have 10 free parameters in the chiral expressions. 
Of these, two appeared previously at tree level: 
the decay constant / and the ratio, fi, of squared 
Goldstone-meson mass to the sum of valence 
quark masses, m x + m y . There are also analytic 
contributions proportional to L^, L$, 2L§ — L5, 
and 2Lq — L4. Finally there are four 0(a 2 ) pa- 
rameters that appear because of taste violation 
[HJ: a 2 5' A and a 2 S' v , which are the taste- violating 
"hairpin" parameters that enter into NLO chi- 
ral logs, and a 2 C and a 2 F, which give analytic, 
0(ma 2 ) contributions. 

Since the data is very precise (~ 0.1-0.4% sta- 
tistical errors) and we include quark masses as 
large as mf ys (to extract K physics), we must 
go beyond NLO to get good fits. We include all 
C(m 3 ) NNLO analytic terms. NNLO chiral logs 
and taste-violating analytic terms are unknown 



and not included. But for larger masses, where 
NNLO terms are non-negligible, the logarithms 
should be changing slowly and thus well repre- 
sented by analytic terms. Similarly, for larger 
masses, 0(m 3 ) terms should be more important 
than taste-violating 0(m 2 a 2 ) or 0(ma 4 ) terms. 
When expressed in "natural" chiral units, the 10 
NNLO 0(m 3 ) parameters (5 for the decay con- 
stant and 5 for the mass) should have coefficients 
of O(l) if chiral perturbation theory (XPT) is well 
behaved. We constrain them to have standard de- 
viation cr = l around using Baysean priors [H]. 

The strength of the chiral log terms is governed 
by the "chiral coupling," l/(16vr 2 / 2 ). Typically, 
one takes f — f, where / is the bare (tree-level) 
decay constant. For better convergence of XPT, 
it seems reasonable to put a physical parameter 
here, e.g., , f = f n or / = (the difference 
is NNLO). In practice, we try 3 approaches: (1) 
/ = /, (2) / = U, and (3) / = U/y/u, where 
lu is a new fit parameter allowed to vary around 
1 with a = 0.1. So far, the best fits have been 
obtained with choice (3); we currently use such 
fits for central values. Choices (1) and (2) are 
now also giving acceptable fits. Only choice (2) 
has been included in systematic error estimates 
below, but choice (1) will be added in the future. 

To test convergence of XPT, a NNNLO term 
in valence quark mass is included for both decay 
constants and masses. Such terms are found to 
have small coefficients; we plan to eliminate them 
for the fits that determine central values. 

We fit both coarse and fine lattice data simul- 
taneously. The 4 taste-violating parameters are 
forced to change by the factor of 0.375 (ratio of 
ct 2 a 2 ) in going from coarse to fine. The remain- 
ing physical parameters are expected to differ be- 
tween coarse and fine lattices by a s a 2 Aq CD ~ 
2%. We therefore also include an additional "scal- 
ing" parameter for each physical parameter. The 
scaling parameter is the fractional difference be- 
tween the physical parameter on the coarse and 
fine lattices. The scaling parameters are con- 
strained to be with a of 0.02 or 0.025 in central 
value fits; this is changed to 0.01 or 0.04 in fits 
used to estimate systematics. 

Finally, the 4 parameters used in smoothing 
n/a are allowed to vary by 1 a. The total is 46 
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Figure 1 . Decay constants with m x = m y . The 
two cyan "fancy squares" are full QCD (m x = 
m y = m q ) points, continuum extrapolated at fixed 
quark mass. The lowest two lines (dotted red 
and cyan) use continuum-extrapolated fit param- 
eters, m x = m y = m q , and either the physical or 
the (fine lattice) nominal strange mass. The first 
confidence level (CL) is computed in the stan- 
dard way; the second treats the Baysean priors 
as if they were additional data points. 

parameters (or 44 with choices (f ) or (2) above), 
although many are tightly constrained. 

Figure ^ shows some decay constants from our 
central-value fit to the partially quenched data 
with m x ,m y <0.8mP hys (194 points). We deter- 
mine physical quark masses from the meson-mass 
part of this fit (see Ref. jl()| 'l with and ttik as 
input; the decay constants are then extrapolated 
to the physical masses and to the continuum. 

Our preliminary results for decay constants are: 

U = 129.3 ± 1.1 ± 3.5 MeV 
f K = 155.0 ±1.8 ±3.7 MeV 
f K /U = 1-201(8)(15) , (1) 

where the first error is statistical; the second, sys- 
tematic. The largest error on and fx is the 
2.2% scale uncertainty. A chiral and continuum 
extrapolation error of w 1.5% has been estimated 



by considering alternative fits. The results agree 
with experiment within errors. 

Preliminary results (in units of 10 -3 ) for the 
Li at chiral scale m, are: 

2L 8 -L 6 = -0.1(l)(tJ) ; L 5 = 1.9(3)(±1) 

2L 6 -L 4 = 0.5(2)(1J) ; U = 0.3(3)(±£) . (2) 

The errors here are dominated by differences over 
The result for 2L§ — L5 is well outside the range 
that allows for m u = |llj . 

Finally, our preliminary results for quark 
masses at scale 2 GeV are mf s — 70(15) MeV 
and m MS = 2.7(6) MeV, where rh is the aver- 
age of the u and d masses. These results use the 
1-loop perturbative result for mass renormaliza- 
tion |12I13| . The rather large error is dominated 
by the 0{a 2 s ) correction to the renormalization 
constant, estimated in Ref. ^2] to be ~20%. 

Additional discussion of the calculation and re- 
sults appears in Ref. [TU| . 
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